TIT/HEP-366/COSMO-83 



Bubble wall perturbations coupled with gravitational waves 



Akihiro Ishibashi f\ and Hideki Ishihara [] 

Department of Physics, Tokyo Institute of Technology, 
Oh-Okayama Meguro-ku, Tokyo 152, Japan 



Abstract 

We study a coupled system of gravitational waves and a domain wall which is the boundary of a vacuum bubble in 
de Sitter spacctime. To treat the system, we use the metric junction formalism of Israel. We show that the dynamical 
degree of the bubble wall is lost and the bubble wall can oscillate only while the gravitational waves go across it. It 
means that the gravitational backreaction on the motion of the bubble wall can not be ignored. 
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I. INTRODUCTION 



Topological defects could be formed when symmetries arc broken in the early universe. Domain walls, the simplest 
topological defect, emerge at the boundaries of true vacuum bubbles nucleating in a false vacuum at first-order phase 
transitions Q. 

Motivated by recent observations which suggest that the present density parameter tto is smaller than unity, a 
possible inflationary scenario which is compatible with an open universe is proposed [Q. The scenario requires two 
inflationary epochs. First, an old type inflation occurs and solves the horizon and the homogeneity problems. It is 
exited through a nucleation of a single vacuum bubble by quantum tunneling. If the tunneling process is represented 
by the Euclidean bounce solution with 0(4) symmetry |J, the interior of the bubble has the 0(3,1) invariance. This 
symmetry naturally induces negatively curved time-slices and an open Friedman-Robertson- Walker (FRW) universe is 
realized inside the bubble. Next, a slow rollover inflation occurs after that. The present value of f^o is determined by 
the form of the inflaton potential. This scenario is called "one-bubble open inflation" . The geometry of the scenario 
is schematically depicted in Fig. 1. 

In the one-bubble open inflation scenario, the observable universe is surrounded by a bubble wall. A conclusive 
verification of the scenario is the detection of existence of the bubble wall. Though the homogeneous isotropic open 
time-slices do not intersect the bubble wall, it is causally visible from the observers inside the bubble. If the bubble 
wall is deformed in some reasons, we may observe the inhomogeneities of the universe produced by the deformation. 
Thus it is especially relevant to investigate the motion of the bubble wall in the one-bubble open inflation context. 

In the thin wall approximation, the motion of a domain wall is described as a three-dimensional timelike hypersurface 
in the four-dimensional spacetime. The small deformation around background domain walls with planar or spherical 
symmetries are represented by a single scalar field on the background wall hypersurface 0]. The equation for the 
scalar field takes the form of the Klein-Gordon equation on the hypersuface. Thus, it is expected intuitively that the 
wave propagation of the scalar field represents the oscillation of the deformed wall. These analysis are based on the 
assumption that all influence of gravity on the wall motion is negligible. 

In the thin wall approximation, the energy-momentum tensor of the wall is concentrated on the three-dimensional 
hypersurface. Thus, the spacetime should be singular at the hypersurface. In the general relativity, the spacetime 
containing such a singular hypersurface should be treated by the metric junction formalism j^] . In this method the 
wall is treated as a source of the gravitational field in which the wall moves. The motions of domain walls with 
spherical symmetry are analyzed by this formalism in Ref . HQ . 

In the metric junction formalism, Kodama, Ishihara, and Fujiwara studied the perturbations of a spherically 
symmetric domain wall system |@J . The interaction between the perturbed wall and gravitational waves is taken into 
account in their work. As a result, it turned out that the deformation of the domain wall occurs only while the 
gravitational wave goes through and does not propagate on the background hypersurface by itself. 

The background configuration of the spacetime considered in Ref. || is constructed in the peculiar way: joining 
two copies of a portion of Minkowski spacetime with spherically symmetric boundary by identifying the boundary at 
the domain wall. So, there is a reflection symmetry with respect to the spherically symmetric domain wall. Thus, it 
is interesting to clarify whether the result is due to the symmetric configuration of the background, or to the domain 
wall's feature itself. This also motivates us to investigate the motion of domain wall on the different background 
configuration: one-bubble open inflation configuration, in which there exists the vacuum energy deference between 
the inside and the outside of the bubble and no reflection symmetry. 

In this paper, we investigate perturbations of a coupled system of gravity and a vacuum bubble in de Sitter spacetime. 
We treat the bubble wall by the thin wall approximation and apply the metric junction formalism to consider the 
gravitational back reaction on the bubble wall motion. We accomplish this in terms of the gauge-invariant perturbation 
method on spherically symmetric background spacetime established by Gerlach and Sengupta J|-[ll]] . Though, in the 
main stream of our study, we consider a bubble wall whose inside is Minkowski and outside de Sitter spacetime as a 
background, we also examine the case that both sides of the bubble are de Sitter spacetimes with different cosmological 
constants. The issue is the extension of that in Ref. 101. 
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In the next section, we specify a background geometry which is used in the one-bubble open inflation scenario by 
using metric junction method. In section 3, we consider the metric perturbations in terms of gauge-invariant variables. 
We find master variables for two modes of gravitational waves and reduce the perturbed Einstein equations into a 
single master equation for them, and give the general solutions of it. In section 4, we write down the perturbed 
junction conditions as the boundary conditions on the master variable. We match the general solution along the 
bubble wall by using the perturbed junction conditions and obtain the global solutions. Section 5 is devoted to 
summary and discussion. 

Throughout this paper, k denotes 8ttG and the unit c = 1 is taken. 

II. BACKGROUND GEOMETRY 

Let us consider a geometrical model in which a spherical vacuum bubble exists in dc Sitter spacetime with a cos- 
mological constant A as a background spacetime. The bubble wall divides the spacetime into two regions: Minkowski 
(M_, cis 2 ,) inside and de Sitter spacetime {M + ,ds 2 + ) outside the bubble. Hereafter, we describe quantities in M_ with 
subscript " - " and in M + , with " + ". 

We assume that the thickness of the bubble wall is small enough comparing to all other dimensions, that is, the thin 
wall approximation. Then the history of the bubble wall is described as a three-dimensional timelike hypersurface 
E. The whole spacetime manifold containing the bubble wall is constructed as M := M_ USU M + by matching the 
boundaries E± := dM± at E. The spacetime metric ds 2 is required to be C° in M. The geometries of the boundaries 
E± are characterized by the induced intrinsic metrics {qijdQ' l dC, J )± and the extrinsic curvatures (Kijd£ % d£ 3 )±, where 
Q. (the indices i, j, ■ ■ ■ run over 1,2,3) are three-dimensional coordinate systems on E±, respectively. The intrinsic 
metrics and the extrinsic curvatures are specified by the junction conditions. 

In more realistic model of the one-bubble open inflationary universe, the interior of the bubble has a slow rollover 
inflationary phase and it should be treated locally as de Sitter spacetime. We will see later that the results for this 
case are essentially the same. For simplicity, we consider the time symmetric configuration as illustrated in Fig. 2, 
instead of a bubble nucleation. 

A. Intrinsic metric and extrinsic curvature 

The background spacetimes have spherical symmetry. In general, the metric with spherical symmetry can be 
expressed in the form 

ds 2 = g^dafdx" = lab {y c )dy a dy b + r 2 {y c )n pq dz p dz\ (2.1) 

where £l pq dz p dz q is the metric on a unit symmetric two-sphere with an angular coordinates z p , i.e., 

il pq dz p dz q = dO 2 + sin 2 9d<p 2 =: dO 2 . (2.2) 

The functions r{y c ) and 7 a b(y c ) are scalar and tensor fields on the two-dimensional orbit space 2 M spanned by the 
two-coordinates y a and each point of 2 M represents a symmetric two-sphere. Hereafter, the indices on the orbit space 
a, b, c, • • • run over 0, 1 and the indices on the unit two-sphere p,q,r,- ■ ■ over 2, 3. 

Let V M , D a , and D p be the covariant derivatives with respect to g^ u , j a b, and Q pg , respectively. The covariant 
derivative is expressed by D a and D p , that is, the Christoffel symbol r£„ associated with the spacetime metric 
g^v is expressed as 

pa pa pr pr 

be DC pq P°' 

r^(^), r; 9 .-r 2 ^(^), (2.3) 
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where T^ c and t r pq are the Christoffel symbols associated with 7 a h and £l pq , respectively. 
The orbit space metrics of M± can take the static form 



where 



(labdy a dy b ) ± = -g{r ± )dtl+g{r ± )- x dr\, (2.4) 



g(r-) = l, g(r + ) = l-H 2 r 2 + , (2.5) 



and H is defined by H := a/A/3. The history of the spherically symmetric bubble wall is described by a timelike 
orbit in 2 M. 

Let (n M )± be the unit normal vectors to the boundaries E±, respectively, such that both {n^)± direct from M_ 
to M + . Since the vectors (n M )± have vanishing z p components, we can identify them with (n a )±. Let (r a )± be 
the future-directed unit timelike tangential vectors to E± on the orbit space. Then the orbit space metric can be 
decomposed as 

(lab)± = -(Ta)±{n)± + (n a )±(n b )±. (2.6) 
The projection tensors (q fiv )± onto £± are represented by 

(q^)± : = (g^)± - {n ll ) ± (n v )± 

= -(t m )±(t v )± + r|0^, (2.7) 

where r M and f^,, are the four-dimensional extensions of r a and £l pq . The four-dimensional extensions of the extrinsic 
curvatures {K pv )± are defined as 

(K^) ± := -(<?XV a O±. (2.8) 



From Eqs. ( |2.7|) and (|2.8|), the extrinsic curvatures (K fjiu )± can be represented by 

{K ab ) ± = -(raTbr^y/nJi, (2.9) 

(K ap )± = 0, (2.10) 

(•Kp,)± = -(rD±r) ± tt pq , (2.11) 

where 

£>/ := r a D a , D ± := u°£> . (2.12) 



B. Background bubble motion 

By using the metric junction formalism ||, the motions of spherically symmetric thin walls are studied in Ref. [0. 
To establish notation, we briefly review the metric junction formalism and determine the motion of a spherical vacuum 
bubble in de Sitter spacetime following Ref. . 

For convenience, we introduce the Gaussian normal coordinate system (x 11 ) — (x>C) m the neighborhood of the 
hypersurface E such that the spacetime metric takes the form 

ds 2 = g^dx^dx" = d X 2 + qijdCdC- (2.13) 

The coordinate axis \ is normal to E and can be set x = on E. 

Let S^v^ 1 ) be the wall stress-energy tensor which is defined by the energy-momentum tensor as 
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Suv(C) '■= um / Tfivdx 
in the thin wall approximation. Then the energy-momentum tensor can be described as 



(2.14) 



(2.15) 



where p — 3H 2 / k is the vacuum energy of de Sitter spacetime and 9(x) is a step function. Provided that the spacetime 
metric g^ u is continuous at x — 0, (<}ij)+ coincides with (<7y)- on S. Hence, 



kij] = o, 



where, for geometrical quantity Q on M±, [Q] denotes 



[Q] := lim{Q| x=+e - Q| X =_J. 



(2.16) 



(2.17) 



With the help of the Gauss-Codazzi equation, the (i, j) components of the Einstein equations are written by 

^R l3 - KK l3 + 2K a K] + d x Kij = k(T 13 - ^T), (2.18) 

where ^Rij is the Ricci curvature on £ with respect to q i3 . Substituting Eq. ( |2.15| ) into this equation ( 2. IS ) and 
integrating it with respect to x m an infinitesimally small interval (— e, +e), we obtain 



Similarly, from the (x, i) and the (x, x) components, we obtain 

{i) D,S 13 = 0. 



where, for Q on M±, Q denotes the mean value 



Q'=\ lin ?,{Q| x =+e + <2lx=-d> 

Z c— >0 



(2.19) 



(2.20) 
(2.21) 



(2.22) 



and ( 3 'Di is the covariant derivative with respect to gy. 

Assuming that a scalar field comprises the bubble wall, we can write the wall stress-energy tensor by a single 
scalar field a as 

= -aq^. (2.23) 
From the condition (2.21), a should be constant. Then, the junction conditions (2.19) and ( |2.20 ) reduce to 

[Kj\ = (2.24) 

and 

K=-?-. (2.25) 
a 



The bubble wall motion is determined by the conditions ( 2.16 ) and ( 2.24 ). Eq. ( |2.25|) is then satisfied automatically. 
The unit tangential vectors to £± are expressed as 



r| :=dy a ± /dr=(i,r)-_ 



(2.26) 



in the static charts (t,r)± of 2 M. Here, r is a propertime along the wall orbit and a dot denotes the derivative with 
respect to r. From the normalization r a r a = —1, it follows that 
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t± =e- 



1 



9(r±) + rl, (e = ±l). 



' g{r±) 

From the orthogonality n a T a = and the normalization n a n a = 1, the components (n , n r )± are written as 



9(r±) ' 



W)± = evs(r±)+r± = 5(r±)i±- 



(2.27) 



(2.28) 



Then, from Eq. (2.11), we get the (p, g) components of the extrinsic curvatures of S± in the form 



\fg{r±) + ^i 



(2.29) 



into the condition (2.24), we obtain 



The (p, g) components of the continuity condition (2. 1£) require r + (r) =r_(r) = r(r). Then, substituting Eq. (2.2£) 

I + r J = ( - ) ■ (2.30) 



where a is the constant defined by 



o 



4kct 



a 



AH 2 + K 2 a 2 ' 



(2.31) 



The value of the constant a gives the minimal radius of the bubble, which is the radius on the instance of the bubble 
nucleation. 



The solution of Eq. (2.3C) is obtained as 



a cosh t. 



(2.32) 



Here and for the rest of this paper, r denotes the normalized propertime by a. Then from Eq. (2.7), the intrinsic 
metric on the wall takes the form 



(2.33) 



q^dCdC = a 2 (~dr 2 + cosh 2 T(2Q 2 ). 

This is the three-dimensional de Sitter spacetime metric with the 0(3,1) symmetry of the bubble wall. 

As a natural extension of the wall metric (2.33), we introduce the hyperbolic charts (t± , \± , @± , 4>± ) defined by 

t_ =: a(x_)sinhr_, y/g(r + ) sinh Ht + =: a(x+)sinliT + , 

r± =■ a(x±)coshr±, (2.34) 

where 



{l/H) sinHx+ for x > — de Sitter side, 
X- for x < — Minkowski side. 



(2.35) 



This is a Gaussian normal coordinate system (2.13) of the considering model. The coordinates x± are related to x as 

X = a(x±)-a, (2-36) 

and hence at the bubble wall 

o(x±)b = «. (2.37) 
In addition, on the wall, it can be set naturally as 

r±=r, 9 ± =9, (f> ± =<t>. (2.38) 
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Then, in the hyperbolic charts, the spacetime metrics take the form 

{ds 2 )± = d X 2 ± +a 2 ( X ±) {-dr 2 + cosh 2 TdQ 2 } . (2.39) 

Since these coordinate systems, which are also called the spherical Rindler charts, do not cover the whole spacetime, 
an analytic extension is necessary to cover the interior of the forward light cone from the center of the bubble, which 
corresponds to the observable universe. By the coordinate transformations, 

X- = «7, r_ = £ - (2.40) 

the metric takes the form 

ds 2 _ = ~dr/ 2 + T) 2 {d£ 2 + sinh 2 Zdfl 2 } . (2.41) 

This is the Milne universe, which has open FRW time-slices. 

In the hyperbolic charts, the extensions of the normal vectors (n tl )± and the tangential vectors (t m )± to x = const, 
hypersurfaces in a neighborhood of S are expressed, respectively, as 

(nn± = (d x± r, (Tn± = ^—AdrV. (2.42) 

a {X±) 

The extrinsic curvatures of the x = const, hypersurfaces are calculated explicitly as 

(Kij)± = -(d x± lna(x±))(fftf)±, (2.43) 

and evaluated on E± as 



1 x/l - H 2 a 2 

{K ij )_ = --(q ij )_, (K lJ ) + = - V — ( qij ) + . (2.44) 



III. PERTURBATIONS 

In this section, we consider the metric perturbations on the background geometry constructed in the previous 
section. We describe the metric perturbations by the tensor harmonics defined in Appendix A. We introduce gauge- 
invariant variables in accordance with Gerlach and Sengupta H and express the perturbed Einstein equations in terms 
of them. Furthermore, we show that the perturbed Einstein equations reduce to a single scalar master equation in 
2 M. 

A. Gauge-invariant perturbation variables 

Let us consider the metric perturbations (h tIil ,dx tJ 'dx' / )± on the background metrics (g tlu dx tJ, dx")±. Due to the 
spherical symmetry of the background geometry, we can expand the perturbations (h fiu dx fl dx u )± by using the tensor 
harmonics on the unit sphere (see Appendix A, for their definitions) as follows: 

for odd modes, 

h$dx»dx" = {r(f {ol) lm)a(V^im) P (dy a dzr + dz*>dy a ) + r 2 (f^ lm )(T^ lrn )„dz*dz«} , (3.1) 

l,m 

for even modes, 

hf v dx»dx v = J2{(f ie) lm)abYi m dy a dy h + r(/< el > lm ) a (V^ i m ) P (dy a dzJ } + dz^dy a ) 

1,171 

+ (r 2 (/ (e V)(T( £C V) pg + r 2 {f^ 2 \ m ){T^ 2 \ m ) pq ) dz^dz*}. (3.2) 
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The subscripts e and o denote the even and the odd modes, respectively. The expansion coefficients, 
f(y c )i fa(y c ), fab(y c ) are scalar, vector, and symmetric tensor fields on the orbit space 2 M, respectively. Hereafter we 
suppress the angular integers I, m and summation J2i m f° r notational simplicity and discuss the modes I > 2, which 
are relevant to the gravitational waves. 

The metric perturbations h flu dx fJ, dx' y have the gauge freedoms. Associated with an infinitesimal coordinate trans- 
formations, a; M — ► x^ + £ M , the gauge transformed metric perturbations are 



h^dx^dx" = h^dx^dx" + Sh^dx^dx" , 

= h^dx^dx" - (V M £„ + V W £ M ) dx^dx". 

The generator ^dx^ of the infinitesimal coordinate transformations is expanded by the tensor harmonics as 

£<>W = r^ o) V( o) dz p , 

^dx" = Z ( a e) Ydy a + r^ e) V^dz p . 



(3.3) 



(3.4) 
(3.5) 



Then, from Eqs. (3.1), ( p.2[ ), and ( |3.3[ ), the gauge transformed expansion coefficients are expressed as follows: 

for odd modes, 



for even modes, 



S f (oi) 

vi? = 

^y(eO) 
f (e2) _ 



-rD 
2 



r 



--&\ 



- D b £ a 
-rD„ ' 



r r 

(e) 



Sp e/i > = --£ (e> . 



(3.6) 
(3.7) 

(3.8) 
(3.9) 

(3.10) 
(3.11) 



According to Gerlach and Sengupta let us introduce the gauge-invariant variables by combining the expansion 
coefficients of the metric perturbations. 



(3.12) 



From Eqs. (|3.6|) and (3/7), we take the odd mode gauge-invariant perturbation variables 

From Eqs. Q-(fIj]), we take the even mode gauge-invariant variables 

T ai := f ab - D a X b - D b X a + i 7afc |/( e °) + 1(1 + l)/( e2 > - ^X c D c r^ , 
where the vector X a is defined by 

X a := r/' el ^ a — -r 2 D a . 



From Eqs. (3.9) and (3.11), we see that the vector X a transforms as 

SX a = 



(3.13) 



(3.14) 



(3.15) 
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B. The perturbed Einstein equations 



The perturbed Einstein equations on de Sitter background spacetime are described by 



- g^V>S*h x ° + H 2 g^h x x - AH 2 h^ = 0, (3.16) 
where □ := V^V^. We express this perturbed Einstein equations by the gauge-invariant variables introduced in 



Eqs. (3.12) and (3.13). Further, we rewrite them in the form of the wave equation for a single scalar field on the orbit 
space. 

For a while, we consider only the de Sitter side (M + , ds 2 ,) and suppress the suffix " + ". All geometrical quantities 
for Minkowski side (.M_,ds?_) are reproduced by taking the limit H — > 0. We shall discuss the odd and the even 
modes separately. 



1. Odd modes 



The odd mode perturbed Einstein equations are written by the odd mode gauge-invariant variables in the form 



A D C ( - r 4 D a I !(/ + -DrT" = 0. 



(3.17) 



D a (rT a ) = 0. 



(3.18) 



From Eq. ( 3.18 ), we can find a scalar variable $( D ) such that 

rT a = e ab D b ^ (o) , (3.19) 

where 

f ab ._ n a T b _ T a n b j g 

the antisymmetric two-tensor. Substituting Eq. ( |3.19| ) into Eq. ( pM7| ), we obtain 

{r 2 D 2 - 2r(D c r)D c - (I - + 2)} $ (o) = 0, (3.20) 
where 2 := D a D a . We can further rewrite this equation in the form 



= 0. 



(3.21) 



This is the odd mode master equation for ^t ) 



2. Even modes 



The even mode perturbed Einstein equations are written by the even mode gauge-invariant variables as the following 
set of equations: 



\{i + ?[% Dc - <^M±^, - (^) (^) ^ - (^) (?f | ^ 

(^) " (^) + D *^' d 



(3.22) 



(3.23) 







Eq. (3.22) implies that there exists a master scalar variable $( e ) such that 

T ab = (D a D b + H 2 lab )<$> [e) . 



Substituting Eq. (|3.24j) into Eq. (|3.23| ), we obtain 

~ u 2 H 



(3.24) 



(3.25) 



For the derivation of Eq. (3.25), see Appendix B. This is the even mode master equation and is the same form with 
the odd mode master equation ( |3.21 ). 



C. General solutions of the master equation 

Let us give general solutions of the odd and the even mode master equations. Here, we omit the subscripts e and 



o, since the both master equations ( 3.21 ) and ( 3.25 ) take the same form. We express the master equations in the 



hyperbolic charts defined by Eq. ( 2.34 ) in the form 

-a(x±?d 2 x± - a( X ±) (d x a( X± )) d x± + d 2 T + 



0. 



(3.26) 



cosh t / \ 1 / ± 

Assuming that (&/r)± can be decomposed into R(x±)T(t), we obtain the equations for R(x±) and T(r), respectively, 

1 



a{x± 

«2 , Ki + 1 



d x± a(x±)d x± + 



14: 



a(x±) 



R(x±) = o, 



coshV^) TW = °< 



where k± are the separation constants. The solutions of Eq. ( |3.27 ) are given by 



a ±ifc±f(x±) 



where 



C(X-) := liix-, 
C(X+) := - In tan 



V 2 



The solutions of Eq. ( 3.2S|) are given by 



T(r) oc F, ±4fc± (-tanhT) 



oc Gi l,l=F«fc± 



1 + tanh r 



where Pi and Gi are the Legendre function and the Jacobi polynomial. 
Let us introduce the set of functions {^i{t 1 k)} defined by 

Mr,k) :=G l (l,lT*fc; 1 + t 2 anhT W fcT 
This satisfies the orthogonality condition 



dr^(r, k')^i(r,k) = 2nd(k-k'), 



(3.27) 
(3.28) 

(3.29) 

(3.30) 
(3.31) 

(3.32) 
(3.33) 



(3.34) 



(3.35) 
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and the set {ipi(r, k)} (—00 < k < 00) is complete in the L 2 (— 00,00) space ol lunctions of r. Then we can express 
the solutions (&/r)± as 

(^j = J dfc^(^^){^±(fc)e ifeC(x±) +B±(fc)e- lfeax±) }. (3.36) 

The mode expansion coefficients A± represent the gravitational waves propagating from the past Rindler horizon 
t = —00 toward the bubble wall and B± the waves propagating from the bubble wall toward the future Rindler 
horizon t = 00, as depicted schematically in Fig. 3. 



IV. GLOBAL SOLUTIONS 

In this section, we first write down the perturbed junction conditions. We can express them in terms of the master 
variables. Then, by matching the solutions of the master equations in M_ and M + along the bubble wall, we construct 
the global solutions of the perturbed Einstein equations in the whole spacetime M. 



The solutions of the perturbed Einstein equations (3.16) should satisfy the following perturbed junction conditions 



[<%] = 0, (4.1) 
1 

2< 



[5K*] = «(ASj " rfSS). (4.2) 



Since the wall stress-energy tensor is given by Eq. ( 2.23 ) with the constant a, the second condition reduces to 



[5K]\ = 0. (4.3) 

As in the case of the metric perturbations, we can expand the perturbed extrinsic curvature SKij by the tensor 
harmonics: 

for odd modes, 

{8K { -% j dx i dx i = r{5K {ol) ) a (y {o) ) p {dy a dz p + dz p dy a ) + r 2 (SK (o2) )(T (o2) ) pq dz p dz q , (4.4) 

for even modes, 

{5K { % 3 dx l dxi = (5K^) ab Ydy a dy b + r{8K {el) ) a {V {e) ) p {dy a dz p + dz p dy a ) 

+ r 2 {{8K^){T^) pq + (5K^){T {e2 \ q } dz p dz«. (4.5) 

For convenience, let us define the decompositions of tensors on 2 M by using r a and n a . A vector V a is decomposed 

as 

V a = -V//T a + V ± n a , (4.6) 

where 

V// := r a V a , V±:=n a V a . (4.7) 
Similarly, a symmetric tensor T a \, is decomposed as 

Tab = T////T a T b + T ±± n a n b - T //± (T a n b + n a r 6 ), (4.8) 

where 

Tim := r a r b T ah , T ±± := n a n b T ab , T /± := r a n b T ab . (4.9) 
We also define the projections of the perturbed extrinsic curvature 
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(SK,)* := g^SK^ 



(4.10) 



which are continuous from Eq. (E73) . 



Now let us express Eq. ( 1.1C ) in terms of the gauge-invariant variables as follows |lQjl] 



for odd modes, 



(5K { f iy ) a 

(«4 o2) ) 



r 



for even modes, 



K e2) ) - 



Dgr 



Tin - -D ± Tl 



Bur 



1 ^ + \{K^ + 2H^X U 



(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 



We note that these quantities contain the gauge dependent variable X±. In the perturbation theory, we treat the 
quantities on the perturbed geometry by embedding to the background geometry. In the present case, we treat the 
geometry constructed by gluing two different spacetimes: Minkowski spacetime and de Sitter spacetime. In addition, 
the geometry contains a singular hypersurface (the bubble wall) at the boundaries of the two spacetimes. We therefore 
must choose the embedding such that it maps the perturbed geometry of Minkowski (de Sitter) spacetime to the 
background Minkowski (de Sitter) spacetime and furthermore the perturbed bubble wall to that on the background. 
Thus, we consider the gauge transformations under the restriction such as £j_ = on the bubble wall. The variable 
X± is gauge-invariant just for the gauge transformations on the bubble wall and can be interpreted as the wall 
displacement variable || . This is discussed in a separate work . 

Using Eqs. (pl|)-(p^), we can rewrite the junction condition (KL3J) in terms of the gauge-invariant variables and 
further in terms of the master variables <E>( ) and ^(e)- We consider the odd and the even modes separately in the 
following subsections. 



A. Odd modes 

First, from Eq. (LI), we see that f^° 2 ^ and /| o1 ^ are continuous. Then, 



From Eq. (3.19), it follows that 



This is rewritten as 



f/, {ol) -\D/,& 2) 



[D±* (o) ] = 0. 



= 0. 



(°) 



1 /*(„) 

2 \ r 



(4.17) 



(4.18) 



(4.19) 
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Second, from Eq. (4.11 



), we see that 



ab n 

e JJh 



r 



-D a rD n - Do $ 



2 *(„) 



Using the master equation (3.21), we rewrite Eq. ( 1.20| ) as 

[*(o)] = 0- 



(4.20) 



(4.21) 



Eqs. ( 4.19 ) and ( 4.21 ) should be imposed on the solutions of the master equation as the boundary conditions at the 
bubble wall. We obtain no new condition from Eq. (4.12). 

Substituting the solutions ( |3.3G| ) into the junction condition (4.21), we get 



/oo 
dki(ji (r, k) {A+ - A- + B+ - B- } = 0, 
- oo 



where 



A± := A ± e ik±c± | s , B± := B ±e - ik±<:± 



E- 



(4.22) 



(4.23) 



Here, e ±lfc ±C± |s are constant phase factors. Multiplying drip^ (r, fc') and using the orthogonality condition ( |3.35 ), 
we obtain 



[A] + [B] = 0. 

Similarly, from the junction condition ( 4.19| ), we obtain 



k — i- 



na<j\ —r 



\a-Us- 



4 I 



)b = o. 



From the conditions (4.24) and (4.25), we get consequently the relations of the mode expansion coefficients, 



B+ = 



-if i 



A 



k 



k + ifx k + ip 



Air, 



(4.24) 



(4.25) 



(4.26) 



where 



(4.27) 



B. Even modes 



From Eq. ( jllj) , the functions /////, f//^ el \ /^ e0 \ and are continuous. Then 

= 0. 



™ = 



(4.28) 



Similarly, from Eq. (4.16), we see that 

[X±] = 0. (4.29) 
It therefore turns out that the vector X a is continuous quantity. Then, from Eqs. ( 2.44| ), (|2.3l| ), and the defini- 



tion ( |3.13| ) of Tab, we get 
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[>>/] = 



///// - (»*/) - 4(K)X ± - \l(l + l)f^ 

-2naX\, 



D/jr 



Hence, 



and equivalently 



It follows from Eq. (4.14) that 



With the help of Eq. (|4.34|) and the relation 



[^} = + 1(1 + l)/( e2 ) + A-J-X/ + 4(K%)X X 

= i[K-]X ± 
= 2naXj_. 



[?////]+ i^)=o, 
[F±±] = 0. 
[T //A ] =0. 



-{KlY + + 2H' = -{Klf_ 1 



Eq. (4.15) yields 



[D±F C C 



(kp) ( r ±± \r t 



= 0. 



By using Eqs. ( |2.44| ), ( |2.31j ), and ( }4.33| ) , we can rewrite Eq. ( KM ) as 



H 



[D ± J*\ = [^ c c ] - kcj - ~?[ 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



Using Eq. (3.22), we see that no new constraint is obtained from the continuity of (bK^) a a . 

Now let us express the junction conditions ( 4.32 ), ( 4.34 ), and ( 4.37 ) in terms of the master variable 3?( e )- It is useful 
to rewrite Eq. ( 3.24 ) in the form 



T nh = rD„D h 



4> 



(e) 



D„ 



D h 



(e) 



(4.38) 



First, by using the even mode master equation ( |3.25 ), we see that the jump of the trace part of Eq. ( 4.38Q becomes 



[T^] = 2 coshr <J y/l - H 2 a 2 ( «9, $(e) 
and the jump of J 7 // ii 



4> 



[J 7 //] = — \ coshr d 2 . + sinhr<9 T 



— { 2 sinh r<9 T 



4>« 



coshr 



(e) 



1(1 + 1) 

2 cosh r 



Substituting Eqs. (4.39) and (4.40) into the condition (4.32), we obtain 



1(1 + 1) 
cosh 2 r 



= -a^ v 7 ! - # 2 < 



4> 



-a 



(e) 



(4.39) 



(4.40) 



(4.41) 
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Next, from Eqs. and ( jOg ), we obtain 



9 T I coshr 



(<0 



= 0. 



(4.42) 



With the help of Eq. (3.25), from Eq. (4.38), we get 



cosh rS 2 . + sinh rd T 



— \ tanhr9 T ( coshr 



1(1 + 1)' 
2 cosh r 



(e) 



G "!)(« + 2) 



0) 



2 cosh r 



0) 



F±± - -Tf. = - { coshra; + sinhr5 r 
2 a 1 



2 cosh t 



4> 



Substituting Eqs. (4.39), (4.42), (4.43), and (4.44) into the condition (4.37), we obtain 



ci- 



te) 



= 0. 



(4.43) 
(4.44) 

(4.45) 



Eqs. ( 4.41 ) and ( 4.45 ) are the boundary conditions at the bubble wall for the even mode master variable 3>( e )- 



Now, let us translate the conditions (4.41) and (4.45) into the conditions on the mode expansion coefficients. 



Substituting the solutions (3.36) into the conditions (1.41) and (4.45), we obtain, for k =^ 0, 



ik ([A] + [B]) = y/l - H 2 a 2 (A+ - B+) + A- - B_, 
A-B = 0. 



Then, using the equation 



1 - v 7 ! - H 2 c 



2/x, 



(4.46) 
(4.47) 

(4.48) 
(4.49) 

i[i k + tfj, 

Eq. ( p~49l ) are the same form as the odd mode case Eq. ( [l.26j ), apart from the difference of the sign of the first 
term. They gives the manner of reflection and transmission of incident gravitational waves by the bubble wall. 



we finally obtain the relations of the even mode expansion coefficients, 

ifi . k 



B 



± 



k 



-A± 



-A, 



One expect that the poles of Eqs. (4.26) and (4.49) might correspond to the characteristic modes of the gravitational 
waves emitted by the bubble wall. Both the poles, in the present cases, are pure imaginary: k = —ifi. The behavior of 
the pole modes should be examined in the Milne universe t_ > r_ in accordance with the one-bubble open inflation 



context. By the analytic extension of <& with respect to the coordinates transformation (2.40), the metric perturbations 
T a b and T a i n the Milne universe are obtained. Then, one can easily observe that T a b and T a corresponding to the 
pole modes diverge at the center r_ = 0. Hence, the solution $ of the pole modes should be excluded from the global 
solutions by the regularity condition at the center. For the same reason, k = even mode should also be excluded. 
As a result, there is no characteristic mode. 

Since, for odd mode perturbations, there is no variable which describes physical deformation of the bubble wall, the 
result for the odd modes is naturally understood. However, for even modes, Xj_ expresses the physical deformation 
of the bubble wall Its behavior is described by the even mode master variables 3>( e )± and hence the bubble wall 



is not a rigid surface. In fact, from Eqs. (4.30) and (4.40), it is obtained that 

1(1 



X, = — 



1 



naa 



cosh rd£ + sinh rd T 



1) 



2 cosh r 



(e) 



(4.50) 



Thus, the result for the even modes is nontrivial. It is concluded that the behavior of the bubble wall is completely 
accompanied with the metric perturbations and the bubble wall has no its own dynamical degree of freedom. 
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It is straight forward to calculate the case that the interior of the bubble is also de Sitter spacetime which corresponds 
to a slow rollover inflation phase in the one-bubble open inflation scenario. In this case we can also obtain the same 
results, Eqs. (4.26) and (4.4S), only apart from the value of a. In this case 

AK(J (4.51) 



16(if? - H 



2 \2 



in 2 a 2 (H 2 + H 2 _) + k 4 (j 4 



where A± := 3H± (ill > H 2 ,) are the cosmological constants of the outside and the inside of the bubble (see Appendix 
C). 



V. CONCLUSIONS AND DISCUSSION 

We have analyzed the coupled system of perturbations of gravity and the bubble nucleating in de Sitter spacetime. 
We have adopted the thin wall approximation to treat the bubble wall. We have solved the perturbed Einstein 
equations on both sides of the bubble and connecting the solutions along the bubble wall by using the metric junction 
formalism. Thus we have obtained the global solutions of the perturbed Einstein equations in the whole spacetime. 
We concluded that the perturbative motion of the bubble wall is completely accompanied with the gravitational 
perturbations; the bubble wall oscillates only while incident gravitational waves go across. The result is essentially 
the same as that in Ref. ||. 

In the context of the one-bubble open inflation scenario, bubble wall perturbations are investigated by several 
authors fl3|| . Recently, Tanaka and Sasaki have shown that, once the gravitational perturbations are taken into 
account, the bubble wall fluctuation modes disappear E3| . Our result is consistent with their work. 

The behavior of the wall clarified in the present work is different from the intuitive expectation such that deformed 
bubble wall oscillates by its own dynamics and the oscillation damps gradually by emission of gravitational waves. As 
shown in Ref. [Q] , the perturbations on the wall can be described by a single scalar field 4> on the hypersurface which 
represents the background motion of the wall: three-dimensional de Sitter spacetime in the present case. When we 
ignore the gravitational back reaction on the motion of the bubble wall, we obtain the perturbed equation of motion 
for the bubble wall as the Klein-Gordon wave equation on the background hypersurface, 



□a 



= 0, 



(5.1) 



where D s is the d'Alembertian on the three-dimensional background hypersurface and the mass term m 2 is expressed 
by H 2 and a. This equation has oscillatory solutions and the non-spherical oscillations of the wall become a source of 
gravitational waves. Thus, if the influence of the gravity might be negligible, the solution of Eq. (|]l]) have described 
the behavior of the wall well. In the present case, however, we see the qualitatively different behavior of the wall 
motion, taking the back reaction of the gravitational perturbations into account. Our results suggest that the back 
reaction of the gravitational perturbations on the motion of a domain wall should not be ignored. 

The difference of our results and usual analysis may appear in the equation of motion for the bubble wall. By using 
the junction condition, we can easily obtain the equation of perturbative motion of the bubble wall coupled with the 
gravitational perturbations as the equation for X±_ on the orbit space. From the junction condition ( [2.2E ), we get 



5K = 0. 



(5.2) 



Then substituting Eqs. (4.13), (4.15), and (4.16) into this condition, we obtain the perturbed equation of motion for 
the wall as 



+ 



- D ,l 
3H 2 



D /f r 



Di, 



1(1 + 1) 



|ff 2 + 3 

I \ KG 

1- 



#2\ 2 q„2„2 
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(5.3) 
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Comparing this equation with Eq. fl5.1[ ), we notice that there are additional terms described by the metric pertur- 
bations. The task which we should do is solving the coupled system of the perturbed Einstein equations and the 
equation of motion for the bubble wall (5.3) simultaneously This work has been accomplished by using the metric 
junction formalism in the present paper. Consequently the quantities of the second line in Eq. (5.2) are comparable 
to the quantities of the first line, i.e., the back reaction of the gravitational waves on the bubble wall motion is not 
negligible. 

The equation of motion (5^3) is obtained from the junction condition which is derived from the Einstein equations. 
It is natural because the Einstein equations contain the equation of motion for matters. We can also derive the 
perturbed equation of motion ( |5.3| ) directly from the Nambu-Goto action for the domain wall. It will appear in 
Ref. @. 

It is interesting to study the quantum fluctuations of the bubble wall in the early universe. Garriga and Vilenkin 
investigated the quantization of the scalar field satisfying Eq. (5.1) fig] . However, if we consider a domain wall coupled 
with gravitational waves, the dynamical freedom of the wall is lost. Thus, another quantization scheme for the coupled 
system may be applied. 

One of the reasons for the result obtained in the present paper may be the fact that the stress-energy tensor is pro- 
portional to the intrinsic metric of the wall hypersurface (see Eq. (2.23)). More general cases are under investigation. 
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APPENDIX A: TENSOR HARMONICS 



Here we give the definitions and the basic properties of the tensor harmonics. We denote the metric of the unit 
two- sphere by 

n pq dz p dz q = d9 2 + sin 2 9d<j) 2 , (Al) 

and the Laplace-Beltrami operator by A2 := QP q D p D q . 

The bases of the tensor harmonics are chosen so that each has a definite parity with respect to the transformation 
(9, <fi) — > (tt — 9, <fi + 7r) which maps each point on the unit sphere to its antipodal point. Under this choice the 
harmonics T which transform as 

T — > (-l)'T (A2) 

are called even, and those which transform as 

T — ► (-l) i+1 T (A3) 

are called odd. 



1. Scalar harmonics 



The scalar harmonics are the spherical harmonic functions Yi m (z p ) which satisfy the equations 

{A 2 + l(l + l)}Y lm =0, (A4) 
d^Yira = imYim. (A5) 
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They are expressed in terms of the Legendre functions as 



w^-J S^M^ rn*.*)*-*. (A6> 



These are all even. 



2. Vector harmonics 

The vector harmonics on the unit sphere are defined by 

(V^i m )p ■= D p Yi m for even modes, (A7) 

{V {o \ m ) p := e pq D"Y lm for odd modes, (A8) 

where e pq is the two-dimensional Levi-Civita antisymmetric tensor. Hereafter, the subscripts (e) and (o) express the 
even and the odd parity, respectively. Any vector on the unit sphere can be expanded by these vector harmonics. By 
definition, the vector harmonics satisfy 

D p (V^ lm y = -l(l + l)Y lm , (A9) 

D P (V (o) lm f = 0, (A10) 

D[ P {V {e) lm ) q] = 0, (All) 

D[ P {V {o) i m ) q] = + l)e pg Y lm , (A12) 

^A 2 + l 2 + I — ij (Vi m ) p = for both the even and the odd modes, (A13) 

The vector harmonics vanish for 1 = 0. 

3. Tensor harmonics 

Any second lank tensor field on the unit sphere can be expanded by the tensor harmonics which are defined by 

{T {e2) i m ) Pq := {AA + + Y lm , (A15) 

{T {o2 \ m ) pq := \(e qr D P + e pr b q )D r Y lm . (A16) 
The harmonic tensor (T^ e0 \ m ) pq is essentially of scalar type and satisfies 

[A 2 + l(l + l)}(T^ lm ) pq =0, (A17) 

(T^ lm )/ = Y lm , (A18) 

D q (T^\ m fP = \{V^ lm y. (A19) 

(A20) 

The harmonics {T^ e2 \ m ) pq and {T^ o2 \ m ) pq are purely tensorial and satisfy 

{ A 2 + I 2 + I - 4} {T^ o2 \ m ) pq = 0, (A21) 

(T (e2 -° 2) lm ) p p = 0, (A22) 

b q (T<> e2 '° 2 \ m yv = ~(l -!)(/ + 2)(V^°\ m y . (A23) 
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APPENDIX B: DERIVATION OF THE EVEN MODE MASTER EQUATION 



We show that the even mode perturbed Einstein equations ( 3.22| ) and ( 3.23| ) reduce to the master equation (3.25). 

(Bl) 



Substituting Eq. ( |3.24| ) into Eq. ( |3.23| ), we get 

(D a D b + H 2 lab ) L = 0, 

where 

L := {-r 2 D 2 + 2r (D c r) D c + (l- + 2)} $ 



(e)- 



(B2) 



Regarding Eq. (Bl) as a second order differential equation for L, we obtain the general solution of Eq. (Bl), denoted 
by L s , as 



L s = dr + VI - H 2 r 2 (C 2 sinhHt + C 3 cosh Ht) , 
where C\, C2, and C3 are arbitrary constants. Thus, <&( e ) should satisfy 

{-r 2 U 2 + 2r (D c r) D c +(l- l)(l + 2)} $ (e) = L s . 
By the transformation of $( e -j such as 



we obtain 



$(e) $ (e) ( / _ 1 )( / + 2 ) + (Z-i)Z(Z +!)(/ + 2)' 



L' := {-r 2 U 2 + 2r (D c r) D c + (l- 1)(Z + 2)} $' (e) = 0. 



(B3) 



(B4) 



(B5) 



(B6) 



The transformation (B5) does not change J- a b- Thus, rewriting $£n by 3?( e ) 5 we can reduce Eq. (B6) to the equation 



for $( e j of the form 



-□2 + 



0. 



(B7) 



APPENDIX C: DE SITTER - DE SITTER CASE 



Let us consider the case that the inside of the bubble is also de Sitter spacetime, where the cosmological constants 

+(-) 



of the outside (inside) are A + (_) = iH 2 ^, , (iJ 2 . < From Eq. ( 2.11 ), the (p,q) components of the extrinsic 



curvatures are 



(K q p )± 



1 - H\r\ 



r± 



~-S q 

u p . 



Substituting Eq. (CI) into the junction condition (2.24), we obtain the equation of the bubble wall motion 



l + f z = 



where 



Ana 



16(#i - H 2 ) 2 + 8n 2 a 2 {H 2 + H 2 ) + K 4 cr 4 



(CI) 



(C2) 



(C3) 



The solution is 
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r = a coshr, 



(C4) 



where r is normalized by a. This conforms to the 0(3, 1) symmetry of the bubble wall. 

As in the Minkowski-de Sitter case, we can consider the metric perturbations of the geometry in terms of the 
gauge- invariant variables introduced in the section 3. In the hyperbolic charts defined by Eq. ( 2.39 ) with 

f 



a(x±) := 7^-siniJd-Xi, a(x±)|s = a, 



we obtain the master equations ( 3.21 ) and ( 3.25 ) for the master scalar $ and solutions 

dkipi{r, k) {A±{k)e lk ^ + B±{k)e- lk ^ } , 



where 



C(X±) := Tin tan 



H± X ± 



The odd mode junction conditions for 5>( \ are described as 



D 



(°) 



1 



(o) 



(o) 



0. 



(C5) 
(C6) 
(C7) 

(C8) 
(C9) 



These are the same forms as the Minkowski-de Sitter case. In terms of the expansion coefficients A± and B± of the 
solution, these conditions reduce to 



B± = 



-A, 



1 -M- 

k + ifi 



k + ip 

For even modes, the junction conditions of the metric perturbation expansion coefficients are 

[^ c c ] = 2naX ± , 

[?////] = - 2k(7X ±, 
[F±x] = 0, 

Hi -H 2 



[T c c ] - KG ( F ±± - -JP 



(CIO) 



(Cll) 
(C12) 
(C13) 
(C14) 

(C15) 



Rewriting these conditions in terms of the even mode master variable <&( e ) an d combining them, we obtain the following 
constraints for 3>( e ), 



,2 , 1(* + 1) 



cosh 



4> 



1 - H 2 a 2 d 



4> 



1 - H 2 _a 2 d 



(l-l)(l + 2) 



These constraints yield the following relations of the mode expansion coefficients 



IjJL 



k + ifj, 



l — 

k + 



(C16) 
(C17) 

(C18) 



Eqs. ( C10| ) and (C18) are the same forms as Eqs. (4.26) and ( 4.49| ), respectively. 
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Open FRW universe 



Slow roll over 
inflatiok 




\ Cosmological 
\ horizon 



Bubble Wall 



False vacuum : old inflation era 




FIG. 1. The Penrose diagram of the one-bubble open inflation universe. 




M + : de Sitter side 



M_ Minkowski side 



FIG. 2. The Penrose diagram of the background geometry. Minkowski spacetime "M_" and de Sitter spacetime "M+" are 
connected at the bubble wall E (thick line). The timelike dashed lines denote \ = const, hypersurfaces and the spacelike dashed 
lines in M_, r\ — const, hypersurfaces, which correspond to open FRW time-slices. 
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FIG. 3. The block enclosed by dashed lines, the Rindler horizons, is covered by the hyperbolic (the spherical Rindler) 
charts. The thick line represents the bubble wall E. The mode expansion coefficient A+(A-) represents the wave propagating 
from the past Rindler horizon in M + (M_) toward the bubble wall E. Similarly, B+(B-) represents the wave from E toward 
the future Rindler horizon in M+(M_). 
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